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I. INTRODUCTION 



The use of digital computers in the real-time control of aircraft has paved the way 
to new aircraft designs that are faster, more maneuverable and safer than ever before. 
With this growing emphasis on digital control design, new tools and techniques have 
been developed to aid in the design and implementation of complex control systems. 
Traditionally, control systems are developed using classical control design methods that 
are costly and time consuming. Newer technologies, like rapid prototyping based on 
modem control methods, integrate the design process and shorten the time required to 
complete a control design from a few years to a few months. 

The purpose of this thesis is twofold: develop the background material for an 
Advanced Control of Aerospace Vehicles course, and design and flight test an altitude 
hold controller for a UAV. 

At the Naval Postgraduate School, students in the Aeronautical and Avionics 
Engineering curriculum receive a full sequence of controls courses with the final course 
being Advanced Control of Aerospace Vehicles, for which this thesis is written. This 
course will introduce the students to the critical aspects of the design, implementation and 
flight testing of the basic controllers for fixed-wing aircraft. The course examines both 
classical (transform) and modem (state-space) control methods for designing complex 
digital control systems. 

The main focus is on the study of sampled-data systems, systems that have both 
discrete and continuous-time components. Devices used for the interface between the 
continuous and discrete components, the analog-to-digital converter and the digital-to- 
analog converter, are covered in detail. This thesis will derive mathematical models for 
each of these operations, and develop tools for analysis and synthesis. 

The class is largely project oriented. Much of the class is concentrated in the 
Avionics/Controls laboratory and the Unmanned Air Vehicle (UAV) laboratory. The 
projects are centered on a Rapid Prototyping System (RPS) developed by the aeronautical 
engineering department for flight testing control algorithms for unmanned air vehicles. 
The system affords a small team the ability to test new concepts in guidance, navigation, 
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and digital control. The RPS consists of commercially available rapid prototyping 
software with an open architecture design to allow for a wide range of applications. The 
application software developed by Integrated Systems Incorporated (ISI), called RealSim, 
allows students to participate in design projects from the initial concept stage to the flight 
testing phase of the design process. This software has two main advantages: 

• The ability to automatically generate higher-language code such as C for the 
designed controller. 

• The system utilizes industry standard I/O devices including digital-to-analog, 
analog-to-digital, pulse width modulation, and serial capability, permitting easy 
connections to hardware. 

The test bed aircraft used is a UAV called FROG , shown in Figure 1.1. This UAV is 
equipped with a complete avionics suite necessary for autonomous flight. 

The accomplishment of this endeavor has led to a number of successful projects, 
including voice controlled flight and an airspeed controller, and is paving the way for 
more complex projects such as autonomous landing. This thesis will describe in detail 
one such design project, the development of an altitude hold controller implemented in 
the FROG using the RPS tools. The automation of the design process is completed in 
three developmental phases: 

• Feedback controller design. In this phase a model of the plant is created. The 
controller is then designed through various methods of classical input/output 
control techniques and/or modem state-space control techniques. The process 
typically involves many iterations to satisfy specific design requirements. 

• Hardware-In-the-Loop Testing. The feedback system is tested with some or 
all of the actual hardware which will be used to control the aircraft. This is the 
final validation of the controller prior to an actual flight. 
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• Implementation and Flight Test. The controller is implemented in the Fight 
Management System used to control the FROG and then flight tested. 

The design process described above is completed entirely with the RealSim series rapid 
prototyping software and real-time control hardware integrated with the system. 




Figure 1.1: UAV FROG 
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II. SAMPLED-DATA SYSTEMS 



The systems and signals studied in the control of aerospace vehicles are referred 
to as sampled-data systems and have both discrete and continuous components, see 
Figure 2.1, where a continuous-time, or analog signal, is to be processed using a 
computer or special-purpose Digital Signal Processing (DSP) chip set. The interface 
between the continuous and discrete systems include the analog-to-digital (A/D) and the 
digital-to-analog (D/A) converters. 
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Figure 2.1: Sampled-Data System 

The conversion of the a discrete to an analog signal by the D/A converter, which 
usually uses a Zero-Order-Hold (ZOH), involves scaling the digital values and converting 
them to a piecewise-constant continuous output. Where the conversion of the analog 
signal to a discrete sequence by the A/D converter is accomplished by a sample-and-hold 
circuit, called a Sampler. 

This chapter will describe the operations of the digital-to-analog converter and 
analog-to-digital converter and methods for obtaining discrete models of continuous-time 
systems. Material presented in this chapter is in the form of classroom notes, for further 
discussion refer to [Ref. 1,2,3]. 



A. DISCRETE MODELS OF SAMPLED-DATA SYSTEMS 

In this section we establish a general method for obtaining the difference 
equations that represent the behavior of the sampled-data system. The resulting 
equivalent discrete system will then be analyzed using the sample-to-sample discrete 
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transfer function of a continuous system between a D/A and an A/D, as shown in Figure 
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2 . 2 . 

Where 



lx = A-x+ Bu 
P = \ 

[y = C -x + D-u 



Pj = 



jx A+ | =Aj-x k +Bj-u k 
In = Cj ■ x k + Dj ■ u k 



Figure 2.2: Equivalent Sampled-Data System 



1. Digital-to-Analog Converters (D/A) 

Digital-to-analog converters usually use zero-order hold or (ZOH): given a 
sequence of samples, u(kT) at t = kT , and holding its output constant at this value until 
the next sample is taken at t = kT + T to produce a continuous signal. The piecewise 
constant output of the D/A is the signal u(t) that is applied to the plant. Thus, the value of 
u(t) consists of steps as seen in Figure 2.3. 
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2. Analog-to Digital Converters (A/D) 

The analog-to-digital converter samples the analog signal y(t) at discrete times and passes 
them to the computer. The time interval between samples y>k is called the sampling 
interval T. The A/D converter has two functions: quantizing the sampled analog signal 
into a discrete set of levels and coding the quantized representation into an acceptable 
format. 
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Figure 2.3: Digital-to-Analog Conversion 



3. Discretization of Continuous-Time Systems 



Consider a linear, continuous-time system 

) x = A - x + B • u 

, 

y - C-x+ D u 

then 



x(t) = e A(l ,o) •*(/„)+ je A(r T) ■ B- u(r)-dt . 



To discretize let: to = kT ; t = kT + T and assuming u is held constant over the sampling 
interval T (ZOH with no delay) 

u(t) = u(kT) ; [kT < t <kT + T\ 



Then 



kT+T 

x(kT + T) = e A(T) • x(kT) + j e A(kr+r - r) ■ B ■ u(r) ■ dr 

kT 
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Let 



then 



£ = (AT + T)-t 
d£ = -dr , 

r = fa => £ = T 
r = kT + T=>% = 0 



Let x(kT+T) = x k+i , then 



** + . = ^' •** + -5-«* -c/£ 



(Eq. 2.1) 



Define 



<ft = e 
r 



AT 



The discrete equivalent of P 



T = \e Ai -B-d£ 



P = { X * +I = °' x * + r ' w * 
17 \y k =C-x k + D-U k 



(Eq. 2.2) 



Note the equivalent discrete system is shift-invariant, since (ft and T are not functions of 
k. Commands to discretize: Matlab— c2d, Xmath—discretize 



EXAMPLE 2.1: Discretize the following CTS 



From equation 2.2 



x = —a • x + a-u 

y = x 



P(s) = 



a 



s + a 



<£> = e 

7 



-uT 



r= Je""* •(a) d£ = (\-e- aT ) 



x M =e-* r -x k +{\ -e-‘ lT )-u k 



y k = x k 
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B. THE Z-TRANSFORM 



1. Definition 

Given a sequence {**}, k = - <x>, . . . oo , let 

oo 

X(z) = ^x k -z~ k , n> < |z| < Ro , (Eq. 2.3) 

*=-oo 

where we assume we can find values of r Q and for which the series converges. 

As in the case of the Laplace transform, Eq. (2.3) is considered an operator that 
transforms a sequence Xk into a function X(z), symbolically represented by 

X(z) = Z{x k } 

The Xk and X(z) are said to form a z-transform pair denoted as 

X(z) Z{x k } 



EXAMPLE 2.2 Consider the sequence given in example problem 2. 1 



a 



P(s) = 

s + a 

By Eq. (2.3) the z-transform of Xk is 



x k = e 



■akT 



k=0,l,2,... 



Z{x, } = X(z) = £ e -'z-* = . z-') 

*=0 A=0 

Note: a infinite geometric sum is given by 

S. = - A- ; H < 1 

;«=o l— a 

Letting a = ( e~ aT • z _l ) 



m = 



e a 1 • z 1 1 < 1 ; \z\ > e 



1 -(e~ aT -z- ] ) 



- aT 



2. Some Properties of the z-Transform 

Listed below are some properties of the z-transform useful to the material 
presented in this section. For reference, a more complete list is found in [Ref. 1]. 
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Linearity: given {x k }, (y k j 



Z(x k + yO = Z(x k ) + Z(y k ) 

Z[a (x k + y/J] = Z(a x/J + Z(a yO = a X(z) + a Y(z) ; where a 



Time Shifting: 

Gives transform pairs 

Demonstrated by 



Z{x(k + n ) } = z- n X(z) 



x(k+l) <-> z X(z) 



x(k-l) <-> z' X(z) 



Z ( x k +i) = z (xO 
Z (Xk+i) = YjX k+ , ■ z~ k = X ** +1 



-I; 



*+i 






= z-X(z) 



Final Value Theorem: 

Limx(k) - Lim(z - 1 )X(z) 

k — >oo z — >1 

if such limits exists. 



EXAMPLE 2.3: Find the z-transform given in example problem 2. 1 . 



Jx k+i =e- uT -x k+ (\-e-“ r Kj 

If* = ** J 

zX(z) = e-‘" ■X(z) + (\-e- aT )-U(z)) Y(z) _ l-e~" 7 
Y(z) = X(z) \ U{z)~ z-e- uT 



constant 
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C. ANALYSIS OF DISCRETE-TIME SYSTEMS 



1. The Discrete Transfer Function 

Consider the continuous-time and discrete-time systems shown in Figure 2.4. 



u(t) 




y(k) 

— ► 



Figure 2.4: System Transfer Functions 



The transfer function for the continuous-time system is 

Y(s) = {c(sl - AT' ■ B+D }• U(s) 
For the discrete-time system compute 

Z J** + i = °** + r ‘ w * j 

It* ~C-x k + D- u k J 



z ( x k + ,) = z (o-jt* +r-«) 

z-X(z) = 0-X(z) + T-U(z) 
x(z) = (zi-d>y'-r-u(z ) 

Z(y*) = Z(C-x k + D-u k ) 

Y(z) = C ■ X(z) + D-U(z) 

Y(z) = {c-(zi-<t>y 1 r + d}-u(z) 

V v / 

1 , 

Discrete transfer function 

^=c-(z/- or'-r+D 

U{z) 

EXAMPLE 2.4: Find the discrete transfer function from example problem 2.1 
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Given: 



P(s) = 



a 



Continuous-Time 



Discrete-Time 



s + a 



T(s) = [C • (si - A)~' -B + D ] • U(s) => T(5) = U(s ) 

s + a 



1 

o = e' aT r = \e' a f • (a)d% = (1 - e ~ aT ) 



From Eq. (2.4) 



U(z) z-e 

Now consider the sampled-data system in Figure 2.5. 



-aT 



- aT 




Figure 2.5: Prototype Sampled-Data System 



u(kT) = 



Let 

fi;* = o 

| 0;£*0 

Then u(k) is the unit impulse. The response to a unit impulse determines the impulse 
response of the D/A converter: 

«(/) = i(/)-l(/-D 

Then 



U(s) = -~ — = -(l-e"') 



s s 



and 



y(t) = L 



J G(s) 



.(,-c-)) 



y(kT) - y(t), (k-\)T<t<kT 



12 



Finally, 



G(z) = Z(y(kl ")) 
r G(s) 



= ZU" 



* 



•(l-e-)l 



= (i-z' 1 )-z^r l 



G(5)' 



= (i-2-')z 

Where the “Z, _/ ” is dropped for convenience. 



V $ 

r G(s)> 

) 



s 



EXAMPLE 2.5: Compute the discrete transfer function for 

a 



G(s) = 



s + a 



Then 



The corresponding time function 



G(s) 1 1 



s s s + a 



y(kT) = \(kT)-e- akr -\(kT) 



The z transform 



4— } = Z(y(kT)) = -i- z — = - 7(1 

[ s J z - 1 z-e (z-l)(z-e ) 



From equation 2.5 



G(z)=(,- z -')l 



(z- l)(z-e") 



1-e 



z-e 



aT 



-aT 



(Eq. 2.5) 
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The transfer function can also be obtained from the system's response to a unit 
pulse. Suppose the system input, iik for k > 0, and the initial conditions are known, then 
writing out each term of the difference equation gives: 



x, = <t> • x 0 + T • u 0 

x 2 = o • x, + r • w, = o • (o • x 0 + r • u 0 ) + r • «, = <t> 2 • x 0 + <t> • r • u 0 + r • », 

x k =<& k -x 0 + 0*" 1 -r-n 0 + o>* -2 -r-w, 



Therefore, 



Suppose xq = 0 and 



then 



In general. 



k - 1 



x k =(I)A -*o +X 0 *’' ' r ' w < 



/=0 



A-l 



y k = CO* x 0 +y, C-O*-' T u, 



/=0 K{k-i) 



fi;A = 0 

jo ;£*0 ’ 



y k = =hj(k), 



/= 0 



y k = £Ay(A-0-«,- 

A=-oo 



(Eq. 2.6) 



2. Block Diagrams and State-Space Descriptions 

The discrete transfer function, in the z-domain, represents a linear algebraic 
relationship, accordingly multiple linear systems may be described by a system of linear 
equations. 
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Parallel Connections: The system response of a parallel combination of two LT1 
systems is the sum of the single-path transfer function, Figure 2.6. 



u 



> H,(z) 



"+ y 

°r~+ 



M) 



Figure 2. 6: Parallel blocks 



Y{z) = [//, (z)+// 2 (z)] • U(z) = H(z) • U{z) 

W) ' 



Cascade Connections: The system response of two LTI systems in series is related by 
the product, Figure 2.7. 



H,(z) 



H 2 (z) 



y 



Figure 2.7: Cascade blocks 



Y\{z)~ //, (z)-U(z)] 
Y(z) = H 2 (z)-Y ] (z)j 



Y(z) = H x {z) ■ H 2 (z) ■ U {z) 



Feedback Connections: The transfer function of a single loop, Figure 2.8, is given by: 



- wrvJLw Hdz) 


y 


J 


k * 

H 2 (z) 





Figure 2. 8: Feedback transfer function 
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